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The  problem  of  scaling  the  phenomena  of  the  collapse  r-f  the 
wake  of  a  self-propelled  body  in  a  stratified  fluid  is  examined.  Conditions 
are  derived  such  that  the  F roude  number  (based  on  body  speed,  length 
and  Vaisala  period)  is  equal  for  both  model  and  prototype.  Conditions 
are  also  obtained  for  the  existence  of  buoyancy  and  inertial  subranges  in 
the  wake  turbulence.  These  conditions  are  applied  to  determine  the 

smallest  size  model  for  which  the  phenomenon  of  wake  growth  and  collapse 
is  properly  scaled  from  the  prototype. 
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In  this  note,  we  consider  a  self-propelled  bod/  of  length  L 
traveling  at  a  speed  U  in  a  stratified  fluid  which  is  characterized  by  a 

******  t  I 

Vaisala  frequency  N  .  The  diameter  of  the  body  is  D  ,  the  diameter  of 
the  propeller  is  ci  and  it  will  be  assumed  that  L/p  ~  10  and  P/d  ~2 
so  that  /- / p  ~  20  .  For  the  full  scale  or  prototype  body,  typical  numerical 

values  are  U  =6  knots  dz  3m/sec,  i.  ■»  / 0  0  tTi  and  M  =  10~^sec"\ 

The  scaling  parameter  is  the  Froude  number 

f  =  U/(NO) 

which  is,  of  course,  equivalent  to  Richardson  -lumber  scaling  since 

Ri  =  1  /  F2 

In  the  usual  way,  jj  ,  D  and  N  are  chosen  so  that 


where  the  subscripts  p  and  M  refer  to  the  values  of  the  parameters  for  the 
prototype  and  model  respectively. 

In  both  the  model  and  prototype,  the  physics  must  be  the  same  for  a 
meaningful  model  experiment.  The  propeller  produces  turbulence.  The  wake 
grows  by  entraining  fluid  with  a  consequent  decrease  in  the  turbulent  energy 
and  increase  in  the  potential  energy  of  the  fluid  within  the  wake.  Simultaneously 
there  is  viscous  dissipation  of  the  turbulence  energy.  Finally,  when  the  turbu¬ 
lent  energy  drops  to  a  sufficiently  low  level,  the  wake  collapses. 

One  of  the  consequences  of  Ko's  model  of  the  growing  and  collapsing 
wake  (Ko,  1971)  is  that  the  time  for  the  wake  to  collapse,  -fc  ,  is  nearly 
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independent  of  the  Froude  number  and  is 

tc  Aa  2/N 

This  is,  of  course,  in  agreement  with  the  available  experimental  data 
(Schooley  and  Stewart,  1963;  Van  der  Watering,  1966). 

Another  time  scale  is  that  of  the  turbulence  (Batchelor,  1967) 

fT-  •*/« 

where  £  is  the  integral  scale  of  the  turbulence  and  i.(  is  the  turbulent 
intensity. 

Clearly  the  ratio  between  these  times  must  be  equal  for  both  model 
and  prototype  if  the  phenomenon  of  wake  growth  and  collapse  is  to  be  cor¬ 
rectly  modeled.  This  ratio  is 

4/fT  =  Zu/(Ni) 

The  turbulent  intensity  is  a  fraction  of  U 

U  =  LI  9  o  <  <  l 

Ko  suggests  that  0,2S  just  behind  the  propeller  and  this  is  in  rough 

agreement  with  the  results  obtained  by  Naudascher  (1965).  It  is  also  to  be 

expected  that  the  integral  scale  is  a  fraction  of  the  propeller  diameter,  d 
i.  e.  , 

l  =  pd  ,  O  <  (3  <  I 

With  these  assumptions 
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and 


/M  . 

.  /  2U 

_  /*4J  \ 

'  [ Ni 

\  Z  \'?Nd  ) 

~4(?1(nd)  "4(p|Ff 


Therefore 


(t),=  (i 

and  only  if 


M 


That  th‘S  —  °f  thC  ttee  f0r  the  »°  COlhp,.  to  the  r.harart.wnH  - 

—  °f  tUrbulc"ce  «h°  same  value  for  both  the  model  and  in- 

Ei^Stype  if  .he  turbulent  intensity  is  the  sa„e  fraction  of  the  body  speed  and 
tte_ip»egral  scale  i.  the  same  fraction  of  the  propeller  diameter  In  hn.h 
and  prototype.  In  both  model  and  prototype  the  self -propulsion  is  accomplished 
by  a  propeller;  both  the  model  and  prototype  propellers  are  geometrically 
similar  and  the  ratio  of  the  rotational  speed  of  the  propeller  to  body  speed  is 
the  same  for  both  model  and  prototype.  Since  both  the  geometry  and  the  velo¬ 
city  triangles  of  the  blades  are  similar,  one  would  expect,  provided  that  the 
model  propeller  produces  a  turbulent  flow,  that  the  ratios  of  the  integral  scale 
to  the  diameter  and  the  turbulent  intensity  ,o  the  body  speed  are  invariant  with 
scale.  That  is,  we  expect  or  and  (3  are  truly  constants.  Then  the  rati, 

(V/T)  is  also  invariant.  Finally,  we  can  estimate  this  ratio.  If  we  take 
=  »/<A  and  p  =  j/2.  , 
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(V-o = 2^yp  -  /*> 

i 

These  estimates  suggest  that  the  collapsing  wake  can  be  properly 
modeled  at  any  scale  provided  that  the  wake  is  turbulent  at  that  scale.  We 
must  therefore  ask,  is  the  wake  turbulent  behind  the  propeller  and,  say,  up 
to  the  point  of  collapse  ?  It  is  clear  that  the  Reynolds  number  will  be  very 
different  in  model  and  prototype,  but  there  must  be  a  similarity  in  that  a 
buoyancy  range  and  an  inertial  range  must  be  present  in  model  scale  if  they 
are  present  in  prototype  scale. 

We  will  therefore  examine  the  flow  at  two  sections;  just  behind  the 
propeller  and  at  the  poi  it  of  collapse.  Taking  ‘he  time  of  collapse  as 

ic  =  2//V 

the  point  of  collapse  is 

=  u-tc  =  zu/n 

Therefore 

(^cM)=  2^r)=  2(z)(aTd)  ~ 120 

using  the  values  of  U  ,  £>  and  H  for  the  prototype.  So  we  must  estimate 

the  state  of  the  flow  at 

/x  =  0  and  4  =  (ZO  d  , 

In  all  of  the  numerical  work  below,  it  will  be  assumed  that  the  speed  of  the 

model,  4.  =  1  m/sec. 
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A  necessary  condition  for  the  existence  of  an  inertial  range  (Phillips, 
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1966)  is  that 


We  will  assume  that,  just  behind  the  propeller  » 


u=  u/4  ,  1  =  J-h- 


With  Lj  "  1  m/sec,  <vL  =  1/20,  and  V  "  10~^cm^/sec, 


with  L,  in  cm.  It  is  clear  that  any  model  longer  than  about  50  cm  easily 

\)z 

satisifies  the  condition  Rto  >>  *  • 

Next  we  will  estimate  at  /X-I20(l  .  We  need  estimates  of  H 

and  Ji  at  ~  12  0  (L  .  We  will  use  the  only  experimental  results  which  we 
have,  those  of  Naudascher  (1965).  Using  the  data  in  Naudascher's  Table  3 

ft 

and  Figure  1 5 


0.5- 


•j  *  (o  3 


*  Naudascher  uses  D  for  the  diameter  of  the  turbulence  generator  and  1— 
for  the  integral  scale;  our  cL  and  j[  ,  respectively. 
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Therefore 


R'L  (AS*<03f(Ud/»  f 

-  (3S  it  lo’3'/*' (U  L />>'/'' (d /l) 


with  L  in  cm.  We  thus  need 


and  it  appears  that  =  100  cm  will  be  adequate.  Therefore  there  will  be 
an  equilibrium  range  at  collapse  if  L  jfc.  1  meter. 

Next,  we  ask  under  what  conditions  a  buoyancy  range  will  exist.  It  is 
clear  that,  just  as  in  considering  Rt  ,  the  crucial  conditions  occur  at 
^^/20d  •  The  buoyancy  wavenumber  Kb  which  separates  the  buoyancy 

range  and  the  equilibrium  range  is  (Lumley,  1964;  Phillips,  1966) 

Kb  *(N3/e  )'A 

where  £  is  the  rate  of  energy  dissipation.  For  a  buoyancy  range  to  exist 

(1/7)  »  Kb 

where 

H  =  (p3/  <=/* 

is  the  Kolmogorov  micro-scale.  Taking 

6  =  u3/.£  , 
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this  condition  is 


uVO^AO  »  I  . 

I  < 

This  can  be  written  as 

(w/u)3  FZ( d/l)(P/cL)( d/l)(u l/„)  »  /  , 

or  with  ll/u  —  7  *  10  •  -0,5  from  Naudascher's  results,  and 

ZJ  =  1  m/ sec. 

L  »  I 

with  L  in  cm.  Again,  it  seems  reasonable  that  /_  as  1  meter  will  be 
adequate.  * 

Finally,  we  might  require  that  the  Reynolds  number  based  on  the 
dissipation  length  scale 

R*  =  (u  \/»  )  »  |  . 

Again  we  consider  only  the  case  "X  ss  120  Ji  .  Naudascher  (Figure  13) 

gives  curves  of  X/A*  (fa/ci  =  0,l5$)  as  a  function  of  V  (v~*o) / U  fa 
Then,  taking  O  , 


~  z/l 


*  At  ~iZb  d  ,  the  ratio  uV £  is  actually  larger  than  £  by  the  amount 
of  energy  extracted  in  buoyancy  range.  Including  this  consideration,  this 
condition  is  given  by 

L  »>  I 

Again,  it  seems  reasonable  that  L  ♦  100  cm  is  adequate. 


with  L  in  cm.  For  /  m 
to  0.  2.  We  will  take 


<  L-  <  10  In  , 


lies  in  the  range  0.  5 


Atid  Xoi/d  =  0,  !3$~ 
to  obtain  an  estimate  of  fix  .  Therefore, 

/?*  ■=  (“/u  )  (  A/A  -  )  ( A-  /d  )  ( d/L.  )(  Li  Li  ») 

-  o.)G>  L 

v/ith  |  in  cm.  Again,  it  appears  bat  L  lm  will  be  satisfactory. 

We  finally  conclude  that: 

a.  (tt  At)  will  scale  correctly  with  the  Froude 
number  p  -U/(N  x>)  if  the  wake  flow  is  turbu¬ 
lent  at  model  scale. 

b.  For  the  model  wake  to  be  similar  to  that  of  the 
prototype,  the  requirements  are: 

for  an  inertial  range,  ^ 

/?/*=  [(u/u)HU)(^/l)  rJ;"»  i  > 

for  a  buoyancy  range, 

(u/u)3m/. f >(pm)(p/u)  i  ; 

and  finally, 

Rx  =  (u/u)(X/Am)(A K/A)(d/L)  »  I 

with 

Ru  =  U  L  /v> 

Using  Naudascher's  results,  we  find  that  in  all  cases  these  inequalities 
are  satisfied  by  taking  =  10^,  i.  e.  ,  U  =  1  m/sec.  l_  =  l  td  with 

V  =  10  cm^  / sec. 
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